Thermodynamics of small superconductors with fixed particle number 
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The Variation After Projection approach is applied for the first time to the pairing hamiltonian 
to describe the thermodynamics of small systems with fixed particle number. The minimization of 
the free energy is made by a direct diagonalization of the entropy. The Variation After Projection 
applied at finite temperature provides a perfect reproduction of the exact canonical properties of 
odd or even systems from very low to high temperature. 
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INTRODUCTION 

Recent progress in single-electron tunneling spec- 
troscopy have revealed the persistence of pairing effect 
even at very small number of particles [H . The tremen- 
dous experimental work in ultra-small metallic grains [ll 
has enabled to systematically investigate the transition 
from large systems, the bulk limit, up to very small sys- 
tems. By varying the number of particles, thermal exci- 
tations or adding external magnetic fields, the smearing 
of superfluid-to-normal phase transition, the survival of 
pairing correlations, the odd-even staggering 0,01 and/or 
possible re-entrant effects Q have been carefully ana- 
lyzed. These studies have underlined the importance of 
finite size effect on pairing correlations and the necessity 
to develop theories beyond the Bardeen-Cooper-Schieffer 
(BCS) or the Hartree-Fock-Bogoliubov (HFB) ones that 
properly account for particle number conservation. Some 
of these studies arc at the crossroad with nuclear physics 
where systems contain very few to several hundreds of 
nucleons Q and some of the approaches that arc used 
nowadays to deal with particle number conservation, like 
projection techniques [7,, ^] have been imported in con- 
densed matter In this case, improvement beyond the 
BCS and/or HFB, is obtained by considering a state with 
good particle I^'at) = P^I^q), where Pjv is the projector 
on N particles while |$o) denotes a quasi-particle (BCS 
or HFB) state. The explicitly breaking of the symme- 
try, the U{1) one in the present case, allows to grasp the 
physics of pairing while its restoration is required to de- 
scribe the onset of pairing in very small systems (see for 
instance Fig. 1 of ref. @)- 

A natural extension of this approach able to provide a 
canonical description of finite system at thermal equilib- 
rium has been proposed already some times ago [lO[ by 
considering a many-body projected density Dpf written 
as (see also [lH): 



or alternative theories have been proposed. In ref. [12 1 



Dn = —PNGXp{~(3h)PN, 



(1) 



where Z = Tr(Pjv exp(-/?/i)PAr), (3 = l/ikeT), and h is 
the quasi-particle effective BCS or HFB hamiltonian. In 
view of the complexity of this approach, approximations 



general projection formalism was developed and largely 
applied in the static-path-approximation. The problem 
of particle number projection at finite temperature was 
also addressed in the context of the thcrmofield dynamic 
[l3j but no applications have been done till now. Starting 
from a mean-field plus pairing description in the Grand- 
Canonical ensemble, several improvements of increasing 
complexity have been proposed to correct from parti- 
cle number explicit non-conservation. Along this line, 
a Modified BCS theory has been introduced where 
part of the statistical fluctuations is directly incorporated 
in the quasi-particle transformation. This approach has 
been further improved by extending the Lipkin-Nogami 
approach to finite temperature, projecting onto good par- 
ticle number after variation or adding quantum fluctua- 
tion associated to RPA modes Note however that 
its justification and applicabilit y es pecially at high tem- 
perature remain to be clarified [l6[. On the other hand, 
starting from a functional integral formulation and treat- 
ing approximately the collective fluctuations around the 
mean-field path, is shown to provide a suitable tool over 
a wide range of t emp eratures but breaks down at very 
low temperature |17| . An approximate scheme to deal 
with quantal fluctuations consists in the use of a Grand- 
Canonical plus a parity-projected technique [3, H, 
which allows to describe qualitatively odd-even effects 
but still suffers from of abrupt and/or spurious phase 
transitions Even in very schematic models jl9| . un- 
less an exact treatment is made either by direct diagonal- 
ization I^O] or by quantum Monte-Carlo techniques [Hj , 
a canonical flnite-T method based on mean-field theory 
and valid at arbitrary small or high temperature remains 
problematic and appears as a challenge in this field 0. 

While the results presented in ref. [lo| were very 
promising, this method has never been applied due to 
its complexity. Here, we apply for the flrst time the 
method proposed in ref. to the Richardson hamilto- 
nian at thermal equilibrium and show that this approach 
provides a proper description of both thermal and quan- 
tal fluctuations from very low to high temperature. The 
canonical description of a quantum finite system can be 



obtained by minimizing the Helmholtz free energy F 

SF = S{Tt[HDn] - TS) = 0, (2) 

where S denotes the entropy associated to the projected 
density ((J), i.e. S = — kB'Ti-{D]\f In Dj^). The approach 
is apphed to the pairing Hamiltonian written as [l9|: 

H = - '^AiBB)cJ„c,^ - G^cJ+cJ_c, _Ci,+ ,(3) 

i,a—± i.j 

where B is an external magnetic field. For not too big 
systems, thermodynamic quantities can be studied in dif- 
ferent statistical ensembles without approximation by di- 
rect diagonalization of the Hamiltonian in different se- 
niority spaces [ioj. 

The results discussed below are obtained for a system 
of = 10 doubly- folded equidistant levels whose energies 
are 

e^=(^t-^in + l)^Ae, t = l,....,n (4) 

and a pairing strength G = QAAe. In the following, 
the total energy, pairing gap and the temperature are 
given in units of Ae. To take advantage of the U{1) 
symmetry breaking, the hamiltonian h is written as a 
sum of quasi-particle excitations h = EkCt\,ak, where 
the Ek denotes the eigenvalues of the underlying HFB 
hamiltonian, while the quasi-particle creation operators 
write 

"I = ""^4,+ ~ ^kCk-, a| = "fe4,- + ^'fcCfe,+ . (5) 

Similarly to what is done in nuclear physics, two lev- 
els of complexity exist in the application of projection 
techniques. The projection can be made either before 
(Variation After Projection [VAPl) or after (Projection 
After Variation [PAV]) variation[6|. The latter is much 
less demanding since it only requires to solve finite tem- 
perature BCS (FT-BCS) equations and make projection 
without minimizing Eq. As an illustration, the tem- 
perature dependence of the energy {E) and the associated 
heat capacity defined through Cy = d{E)/dT obtained 
with FT-BCS (dashed line) and FT-PAV (dotted line) 
are compared to the exact result (thick line) in figure 
[T] for 10 particles. The exact solution is obtained 

following ref. [1^. As it is well know, in addition to 
the systematic overestimation of the energy, the FT-BCS 
theory suffers from the sharp superfluid to normal phase 
transition as the temperature increases. On opposite, the 
exact solution display a much smoother behaviour. It is 
clearly seen in this figure that, except in the very small 
temperature case, the FT-BCS4-PAV does even a worse 
job and does not cure the threshold effect. 

Extrapolating the improvement generally observed at 
T = (9[ to the finite temperature case, one can antici- 
pate a much better description if VAP is performed. In 
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FIG. 1: Evolution of the energy (Top) and heat capacity (Bot- 
tom) obtained with the FT-BCS (dashed line), PAV (filled 
circles) and exact solution (thick line) for a system of A'^ = 10 
particles. 



that case, the variational principle ^ should be mini- 
mized by both varying the components {uk,Vk) and the 
energy E^ consistently [lo| . While in principle possi- 
ble, such minimization has never been performed due 
to the fact that the hamiltonian h and the operator 
Pn exp{—(3h)PN do not commute and therefore cannot 
be diagonalized simultaneously. As a consequence, while 
a guideline of practical implementation has been pro- 
posed long ago in [l3] , except in the case of the two- level 
degenerate system, the predictive power of VAP at finite 
temperature (called hereafter FT- VAP) has never been 
attested. 

In the present work, we applied the FT- VAP following 
the strategy proposed in ref. [loj . In practice, the varia- 
tional principle is minimized by writing first the energy in 
terms of the one- and two-body density of the projected 
density, both of them written as a non-trivial function 
of the Ufc, Vk and Ek (see Eq. (36) in ref. [l3l)- The 
minimization is carried out via a sequential quadratic 
programming method by using the Vk and Ek as varia- 
tional parameters. To compute the free energy without 
approximation, at each iteration of the minimization, the 



entropy is calculated by 



N 



(6) 



where Df are the eigenvalues of the statistical operator 
exp{—j3h)/Z in the Fock space composed by all the many- 
body configurations with N particles. Each configura- 
tion is characterized by 77 pairs and / unpaired particles, 
with 27] + 1 = N. Moreover, since states with a different 
number of unpaired particles cannot be connected by the 
operator exp(— the problem is reduced to the diago- 
nalization of block matrices for each allowed seniority I. 
The required computational cost is thus given essentially 
by two operations, i.e. the calculation of the matrix ele- 
ments of the statistical operators and the diagonalization 
itself. For the latter, a standard QR algorithm is used. 
The calculation of the matrix elements is done by using 
the bit representation of the many-body states (see for 
example ref. 22|). Each configuration is identified by an 



integer word whose bits correspond to the single particle 
levels and have value 1 or depending on whether the 
level is occupied or empty. In such a way all the matrix 
elements can be obtained by using very simple logical op- 
erations which allow to perform calculations much faster. 

In figure m the result obtained in FT-VAP is compared 
to the exact solution for a system of = 10 particles at 
various temperature. In FT-VAP, the gap is given by Eq. 
(42) of ref. ^ while in the exact case it is computed 
through: 



A = V-G{E - Eo) 
where E is the total exact energy and Eq is given by 

G 



(7) 



(8) 



containing both the single-particle and the self-energy 
terms, rii being the occupation numbers. In this figure, 
we see that, except for the small systematic difference 
observed for the gap, the FT-VAP approach provides a 
perfect description of the thermodynamics of a system 
with fixed particle number in any range of temperature. 
None of the limitations [isl . Kt\ appearing in other mean- 
field based theories are seen. In particular, the entropy, 
that is an approximation in FT-VAP, perfectly matches 
the exact one. The same quality of agreement is found 
also at higher temperature (up to T= 10). 

We further investigated the applicability of FT-VAP 
for odd number of particles. Taking advantage of the fact 
that the FT-BCS density mixes up odd and even parities 
as soon as a non-zero temperature is applied, we used the 
same technique as in the even case. The only difference is 
that now the projector entering in the density (Eq. ([T])) 
corresponds to an odd number of particles. In top panel 
of figure [21 the pairing gap obtained in FT-VAP for N = 
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FIG. 2: Illustration of the predictive power of the FT-VAP. 
From top to bottom, the energy (H), the entropy, the heat ca- 
pacity and the average gap obtained witli the FT-VAP (filled 
circles) are compared to the exact (thick line) and FT-BCS 
case (dashed line) for a system of A'^ = 10 particles. 



10 and A^ = 11 particles is compared to the exact case. 
In bottom panel of this figure, the spin susceptibility x 
defined as la 



xiT) 



52 \nZ 



(9) 



B=0 



is shown for the two cases. 

In the limit of small magnetic field, the susceptibility 
identifies with the fluctuation of the magnetization M — 



'^^BY.a.^^^laCi,a |2l| , i.e. 



X{T) = -^[{M^)-{M? 



(10) 



In small systems, large differences are observed in the 
thermodynamics of odd and even systems [Ij. This is 
clearly seen especially at low temperature for the gap. 
The spin susceptibility further underlines the differences. 
The FT-VAP perfectly grasps the thermodynamics of 
odd systems and one cannot distinguish its result from 
the exact solution. 

In the present letter, we applied for the first time the 
variation after projection approach to describe the canon- 
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FIG. 3: Evolution of the mean gap (upper panel ) and spin 
susceptibility (lower panel) for A'' = 10 (filled circles) and N = 
11 (open squares) particles as a function of T obtained with 
the FT-VAP. The corresponding exact result for the even and 
odd systems are presented respectively with solid and dashed 
lines. In this figure, the spin susceptibility is normalized by 
the bulk high temperature value Xb = 2fj,%/Ae. 



ical properties of a superconducting system. The min- 
imization of the free energy is made with no approxi- 
mation on the entropy. The FT-VAP provides a per- 
fect reproduction of the exact result in the Richardson 
Hamiltonian case both in the low and high temperature 
limit and does not have the limitation of other mean-field 
based approaches. Due to the necessity to make use of 
explicit diagonalization for the entropy, the present ap- 
proach is still restricted to rather small number of par- 
ticles. Nevertheless, we believe that the result obtained 
here is sufficiently promising that in the near future, an 
effort should be made to promote the FT-VAP and make 
it more versatile. It should be mentioned that the present 
method provides a natural extension of the FT-BCS or 
FT-HFB theory presently used to describe nuclei within 
the Energy Density Functional framework applied at fi- 
nite temperature [23|. 
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